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Abstract 



The three flavor linear sigma model is studied as a "toy model" for un- 
derstanding the role of possible light scalar mesons in the tttt, ttK and ttt] 
scattering channels. The approach involves computing the tree level partial 
wave amplitude for each channel and unitarizing by a simple K-matrix pre- 
scription which does not introduce any new parameters. If the renormalizable 
version of the model is used there is only one free parameter. While this highly 
constrained version has the right general structure to explain tttt scatteirng, 
it is "not quite" right. A reasonable fit can be made if the renormalizability 
(for the effective Lagrangian) is relaxed while chiral symmetry is maintained. 
The occurence of a Ramsauer Townsend mechanism for the /o(980) region 
naturally emerges. The effect of unitarization is very important and leads to 
"physical" masses for the scalar nonet all less than about 1 GeV. The ao(1450) 
and -K"o(1430) appear to be "outsiders" in this picture and to require addi- 
tional fields. Comparison is made with a scattering treatment using a more 
general non-linear sigma model approach. In addition some speculative re- 
marks and a highly simplified larger toy model are devoted to the question of 
the quark substructure of the light scalar mesons. 



PACS numbcr(s): 13.75.Lb, 11.15.Pg, ll.80.Et, 12.39.Fe 
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I. INTRODUCTION 



In the last few years there has been a revival of interest [|TJ- |27j in the possiblity that 
light scalar mesons like the sigma and kappa exist. This is a very important but highly con- 
troversial subject. The difficulty is that one must demonstrate their existence by comparing 
with experiment, believable theoretical amplitudes containing the light scalars. However, 
the energy range of interest is too low for the systematic perturbative QCD expansion and 
too high for the systematic chiral perturbation theory expansion |[28|| . Clearly, chiral sym- 
metry should hold but it seems unavoidable to fall back on model dependent approaches. 
Qualitatively, the dominance of tree amplitudes is suggested by the jj- expansion and 
it has been shown by the Syracuse group ||[12],[14],|24| that this approach can be used to 
economically fit the data in the framework of a non-linear chiral Lagrangian which includes 
vectors and scalars in addition to the pseudoscalars. Many related approaches have been 
discussed by other workers ]50 . To put the problem in historical perspective, the theoretical 
treatment of meson-meson scattering has been a topic of great interest for about forty years 
and has given rise, among other things, to chiral perturbation theory and string theory. Nev- 
ertheless, the problem itself of explaining light meson scattering amplitudes from threshold 
to (say) about the 1.5 GeV region is still not definitively solved. Of course, if the existence 
of light scalars is true, it will be a crucial step forward. 

In such a situation, it is often useful to increase one's perspective by studying simplified 
"toy models" . The classic chiral symmetric model which contains a scalar meson is the Gell- 



Mann Levy two flavor linear sigma model [51]. At tree level it yields essentially the same 
7T7T scattering length which is the initial approximation in the chiral perturbation scheme. 
However, compared to that scheme, which uses a non-linear Lagrangian of pions only |32 



it is less convenient to systematically implement corrections. Nevertheless it does contain a 
light scalar meson and it does provide the standard intuitive picture of spontaneous chiral 



symmetry breaking. Furthermore, it is likely to be [55J an exact model close to the QCD 
chiral phase transition. Of course there is an enormous literature on the application of 
the two flavor linear sigma model to tttt scattering. Recently, Achasov and Shestakov Q 
have shown that a qualitatively reasonable picture emerges at the lower part of our energy 
range of interest |34|] by using a scheme which is equivalent to what we may call "K-matrix 
unitarization" . Namely, in the standard parameterization ]5S] of a given partial wave S- 
matrix: 



' = £§-i + *r. (i.i) 



we identify 
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K = T tree . (1.2) 

T tree is the given partial wave T-matrix computed at tree level and is purely real. Such 
a scheme gives exact unitarity for T but violates the crossing symmetry which T tme itself 
obeys. 

For a more realistic application to mi scattering (i.e. inclusion of the /q(980)) as well as 
to tiK, 7T7/ scatterings etc. it is highly desirable to extend this calculation to the three flavor 
case. That is the purpose of this paper. We will see that it provides a very predictive and 
reasonably successful model which gives interesting new insights. 

The three flavor linear sigma model [3(| is constructed from the 3x3 matrix field 



M = S + i(p (1.3) 

where S = S* represents a scalar nonet and = a pseudoscalar nonet. Under a chiral 
transformation q L — > U^q^^ qn — > U^r of the fundamental left and right handed light quark 
fields, M is defined to transform as 

M — ► U L MU R . (1.4) 



To start with, one may consider a general non-renormalizable |37| Lagrangian of the form 



£ = ~Tr (d^d^) - ^Tr (d^Sd^S) -V - V SB , (1.5) 
where Vo is an arbitrary function of the independent SU(3)l x SU(3)r x U(1)v invariants 
h = Tr (MM f ) , I 2 = Tr (MM ] MM ] ) , J 3 = Tr ((MM 1 ) 3 ) , J 4 = 6 (detM + detM f ) . 

(1.6) 

Of these, only J 4 is not invariant under U(1)a- The symmetry breaker V$b has the minimal 
form 

V SB = -2 (A 1 S 1 1 + A 2 S* + A 3 Si) , (1.7) 

where the A a are real numbers which turn out to be proportional to the three light ( "current" 
type) quark masses. 

The Lagrangian Eq. (|1.5|) contains the most general "potential" term Vq but still has 
the minimal "kinetic" term. It is possible [] to also include non-renormalizable kinetic-type 



"See for example, section IV of [38] 
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terms like Tr (d^Md^MM^, Tr (d^MM^d^MM^ + h.c, etc. We shall disregard such 
terms in the present paper. It is interesting to note [37] that the results of "current algebra" 
can be derived from Eq. ( |1.5| ) without knowing details of Vq, just from chiral symmetry and 
the assumption that the minimum of V = Vq + Vsb is non-zero; specifically the "vacuum 
values" satisfy 



[SI) = "X (1-8) 
The "one-point" vertices (pseudoscalar decay constants) are related to these parameters 

by 

F w = ai + a 2 , F K = a x + a 3 . (1.9) 
In the isotopic spin invariant limit one has, 

Ai = A 2 , a± = a 2 (isospin limit). (1-10) 

Many, though not all, of the "two-point" vertices (particle squared masses) may be 



calculated by [37] single differentiation of two "generating functions" which express the 



chiral symmetry of Vq and also using 



\9Sl 



0. (1.11) 



For example, one finds 



m 2 



(n+)=2^^, m 2 (K+)=2^^. (1.12) 



The formula for the mass of the rj' (and of the particles r] and 7r° with which it may mix) 
also involves the quantity 

Many of the three point and four point vertices may be obtained by respectively two times 
and three times differentiating the above mentioned generating equations. The specific terms 
needed for our subsequent discussion are given in Appendix A. 

The present model requires us (in the limit of isospin invariance) to specify the five 
parameters, A ly A 3 , a\, a 3 and V4. These may be obtained by using the five experimental 
input values: 
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m w = 0.137GeV , m K = 0.495GeV, 
m, = 0.457GeV , m n > = 0.958GeV, 

F 7r = 0.131GeV (Inputs). (1.14) 

This is a reasonable, but clearly not unique choice for the inputs. (For example, Fk might 
be used instead of m v ). With these input parameters there are two immediate predictions 
for pseudoscalar properties: 

F K 

P = 2.O5°, -^ = 1.39, (1.15) 
F n 



where the pseudoscalar mixing angle, 9 p , is here defined by 

L — sin£L \ / rio \ 

(1.16) 





cos# p — sin^p 
sin^p cos^p 

r) and rf are the "physical" states while the "unmixed" states are rjs = (<f>\ + <\>\ — 20jj)/\/6 
and r] = (<f)\ + 4>2 + 03 ) / v3- The predictions in Eq. (|1.15|) are qualitatively reasonable but 
not very accurate; the usually accepted value for 9 P , while small, is around —18° and the 
experimental value for ^ is about 1.22. (It is likely that the inclusion of non-renormalizable 
kinetic terms mentioned above will improve this aspect). 

Now, the scalar meson masses are of the most present interest. Analogously to the 
pseudoscalars (tt, K, r/ , r] s ) we denote the scalars by (a , k, cr , cr 8 ). The predicted squared 
masses from the "toy" Lagrangian in Eq. (|1.5| ) are: 

(a 2 - ai) m^ ARE (a+) = 2{A 2 - A ± ), 
2 , + , _ 2(^3 - A x ) 

m BAREl K ) ■ 



«3 — «i 

-4are( K °) = 2iA3 - A >\ (1.17) 
a 3 — a 2 



In this case, the masses of <To, (Jg and their mixing angle 6 S [defined analogously to Eq. Q1.16D 1 
are not predicted. In the isotopic spin invariant limit, which we shall adopt here, the a$ mass 
is not predicted (although it may be reasonably estimated |4(| by taking isospin violation into 
account). Note that we have, in contrast to the pseudoscalar case, put a subscript "BARE" 
on each scalar mass. This is because the pole positions in the pseudoscalar-pseudoscalar 
scattering amplitudes corresponding to scalar mesons may be non-trivially shifted by the 
unitarization procedure of Eqs. (|1 . 1| ) and (|1.2j) . We consider the unitarization to be an 
approximation to including all higher order corrections. Then, in the usual field theoretic 
way of thinking, the pole position determines the physical mass, while the tree level h-bare 
has no clear physical meaning. 
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The tree level irn scattering amplitude is easily computed |37j from Eq. ( [L.5D in the 
present scheme. It involves a four point "contact" amplitude and a and a' exchange dia- 
grams. The resulting form [j] turns out to be remarkably simple: 

,2\2 



A (s, t, u) 




m 2 + cos 2 '?/' 



1™BARE \?) 



m 



BARE 



BARE 



>') 



KareK) - s 



~ ^bare(^) 



- m BARE W 

}. 



;i.i8) 



where s, t and w are the usual Mandelstam variables. The angle ip is defined by the trans- 
formation: 




cosijj —sin?/' 
simp cosijj 

where Q ip is related to the angle 9 S [defined analogously to Eq. Q1.16Q 1 by 




1.19) 



cosV' 



1 

7^ 



(cos6L — v^sin^. 



;i.2o) 



which translates to ip ~ 9 S + 54.7°. With the Lagrangian Eq. ( |1.5|) the amplitude Eq. (|1.18 
depends on the three unknown parameters m B ARE(c), ^bare(c" / ) and ip. 



We can increase the predictivity of the model by restricting the potential V$ in Eq. ( |1.5|) 
to contain only renormalizable terms. The resulting model is the one usually considered 
since it allows for a consistent perturbation treatment (although the coupling constants are 
very large). In any event, we will be working at tree level and "simulating" higher order 
corrections by the K-matrix unitarization procedure. Note that all the formulas gotten above 
with general Vo continue to hold in the renormalizable model; there will just be additional 
restrictions. The renormalizable potential may be written pl[ as: 



V (renormalizable) = [Vt - V u (s a a 2 )] h + -lii(/i) 2 + V 2 I 2 + V 4 I 4 , 
where we have used the notation 



V„ 



'dVo 

dl n 



ab 



' d"V 
dLdh 



1.21) 



;i.22) 



^The sign of A(s,t,u) is the negative of the one in the convention of [37] but in agreement with 
those in fIMl|23l. 



^Note that neither ifi nor 6 S are defined in the same way as 6 S in Eq. (3.6) of [II 
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FIG. 1. Dependence (left) of wj-bare^') in GeV and (right) of the octet-singlet 
mixing angle 9 S in degrees on ^bare(c") ; i n the renormalizable linear sigma 
model. 



As discussed | in |[41|| , we may determine V\ and V2 from the extremum equation Eq. (|1 . 1 1|) 
while V\\ may be expressed in terms of ttt-bare^ ) • Thus specifying m B ARE(o") determines 
the model parameters completely. Actually, tbbare^o) is fixed to be 0.913 GeV just by 
requiring renormalizability, independent of the choice of tbbare^)- Using Eq. (pL - 17]) we 
find that w^bare(^) — 0.909 GeV, independent of whether we make the renormalizability 
restriction or not. Finally, the dependences of h-bare( < ^') and 6 S on the choice of mBARE(c) 
are displayed in Fig. [I]. Choosing the convention where msARE^) < ^bare(0"')> the model 
does not allow for ttt-bare^ ) greater than about 0.813 GeV. Furthermore weare^') must 
be greater than about 0.949 GeV in the renormalizable model. 

In section II we study the simpler two flavor model both for the purpose of review and 
for introducing our notation and the method we will use in the three flavor case. We will 
also illustrate just how well the amplitude can be approximated by a pole in the complex 
s-plane plus a constant. A number of new remarks are made. Section III contains a detailed 
discussion of the s-wave tttt, ttK and 7rr] scattering amplitudes in the unitarized three flavor 
model. Both plots of the predicted amplitudes compared with experiment and numerical 



^Please notice the relevant typographical errors in |41]: 1) In the first of Eqs. (2.2), A3/1V — 1 
should be replaced by A3/1V — A\. 2) In the second line of Eq. (2.5c) AawVu should be replaced 
by Aa 2 wV\\. 3) In the numerator of Eq. (2.8) the factor (4a)~ 2 should be replaced by (2o)~ 2 . 
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calculation of the pole parameters will be seen to be useful for understanding the dynamics. 
A summary and discussion of the calculations of the scalar meson parameters are presented 
in section IV. Section V contains a more speculative discussion on the question of the "quark 
substructure" of the light scalars. It is pointed out that there is a difference in describing 
this at the "current" and "constituent" quark levels. Also, while the linear sigma model is 
set up on the "current" quark basis, it does not uniquely describe the quark substructure. 
In the present model, the initial "current-quark" meson field leads to constituent type states 
which are modified both by details of symmetry breaking and by unitarization. The possible 
richness of the scalar meson system for further study is illustrated by the introduction of a 
larger toy model which includes two different M matrices. 



II. TWO FLAVOR LINEAR SIGMA MODEL 



It seems useful to first review the two flavor case and to make some additional comments. 
We start by exploring the difficulty with a conventional extension of the tree level amplitude 
beyond the threshold region. This also provides the usual motivation for the introduction 
of the non-linear sigma model. 



A. Standard unitarization procedure and its problems 



It is easy to get the two flavor tyty scattering amplitude by taking a suitable limit of the 
three flavor amplitude given in Eq. ( |1.18| ). We simply decouple the a' by setting the a — a' 
mixing angle ip to zero, as is evident in Eq. ( |1.19[ ). Then a becomes Sl ^ 2 , while a' = Sf 
does not belong to the SU(2) theory and decouples; we are left with the tree amplitude 

™BARE \<n ~ S 

The pole term in the second bracket represents the a exchange Feynman diagram. Naively 
one would expect this term by itself to describe a dominance of the low energy amplitude. 
However, the (—1) piece, which comes from the four point contact interaction, is needed 
in this model to satisfy chiral symmetry. It is easy to see that there is a dramatic partial 



A{s,t,u) = —j 



m BARE W - m n 



**Since this formula was gotten as a limit of the SU(3) model with an arbitrary (not neccessarily a 
fourth order polynomial) potential, we see that the tree result Eq. (2.1) is independent of whether 
the SU(2) linear sigma model potential has only renormalizable terms. 
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cancellation of the two terms near threshold. For example, if we take the single unknown 
parameter in the model, m B A RE (a) to be 1 GeV, then the pole term at threshold [s t h = 4m^] 
is about 1.06. This gets reduced to just 6% of its value after adding the constant piece. Near 
threshold we can approximate ^bare(°") ^ [ m ii s ] to get the famous "current algebra" |42| 
formula 

2 (s — ml) 

A(s,t,u)a \ 2 n> - (2.2) 

7T 

We have just seen that this is a small quantity which has arisen from a partial cancellation of 
two relatively large terms. Now if we wish to use Eq. (|2.1[) away from threshold we run into 
the problem of an infinity arising when s = rn BARE {cr) . A standard unitarization procedure 
to avoid this problem would correspond to making the replacement 

> (2.3) 

^bare(^) - 8 ™,l AKE (a) -s- ^bare^) ' 

where T is a width factor. The trouble is that the delicate partial cancellation with the 
contact term is now spoiled near threshold and consequently there will be a very poor 
agreement with experiment in the threshold region. 

The most popular alternative treatment introduces a non-linearly transforming pion field 
and no a at all. (Formally it may be gotten by "integrating out" the a of the linear model 
but this is not the most general formulation). Then the current algebra formula Eq. Q2.2| ) 
is obtained directly from a derivative type four point contact term (as opposed to the non- 
derivative type in the linear model). This approach forms the basis of the chiral perturbation 
scheme (of pions only). The next order correction will involve more powers of derivatives 
and hence will not drastically modify the already reasonable current algebra result. 

A sigma-type particle can be introduced in a general way (independent of the linear sigma 
model) in the non-linear framework by using a standard technique P3]| . In this approach 
the aim couplings are inevitably of derivative type so the a-pole contribution is small near 
threshold and does not drastically alter the current algebra result. This is clearly convenient 
since a regularization of the type Eq. ( |2.3|) will not now alter the threshold behavior dras- 
tically. However, this does not neccessarily guarantee good experimental agreement away 
from threshold. 

It seems worthwhile to emphasize that both the linear and non-linear sigma models 
represent the same physics - spontaneous breakdown of chiral symmetry. The choice of 
which to use is hence primarily a question of convenience in extending the description away 
from threshold. In this paper we focus on studying the linear model, regarding it as a "toy 
model" useful for increasing our understanding. 



9 



To go further, we need the partial wave projection ot the amplitude Eq. ( |2.1|) . Here we 
specialize to the 1=0 projection: 



A 1=0 = 3A (s, t, u) + A (u, t, s) + A (t, s, u) . (2.4) 

The angular momentum / partial wave elastic scattering amplitude for isospin / is 

1 r 1 

T l I {s) = -p{s) dcostf P ; (cosfl) A\s,t,u), (2.5) 

2 7-i 

where A ! (s, t, u) is the isospin I invariant amplitude, 9 is the center of mass scattering angle 
and 

p{s) = (2.6) 

1071^3 

with q(s) the center of mass momentum for, in general, a channel containing particles a\ 
and a 2 : 

T/ is related to the partial wave S-matrix, S l j by Eq. ( |1 . 1| ) . For understanding the properties 
of the a-meson, the Tq amplitude is clearly the most relevant. Using Eq. (|2.1| ) and Eqs. 
( p.4|) - (p.7|) we get the tree approximation 

TLee(s)=a(s)+ (2.8) 

where 

2 /™2 



a (s) = p(s 



™<l AKE {cr) -m 2 



_ 1Q + 1 ^BARE(° r ) ~ m?r ln [ m BARE( ff ) + S ~ 4n V 
S - 477V 2 V ^AReM 



^) = ^#Kare(^)-^ 2 ) 2 . (2.9) 

Note that a(s) in Eq. ( |2.9| ) does not blow up when q 2 = s ~ 4 ^ — > 0. 

Using the partial wave amplitude Eqs. ( |2.8| ) and ( |2.9| ) it is straightforward to give a 
more detailed discussion of the difficulty of regulating the infinity at s — ^bare! ") while 
still maintaining the good agreement near threshold. Consider replacing the denominator in 
Eq. (|2.8| ) according to the prescription Eq. ( p.3|) . The effect of different constant widths V 
in Eq. ( j2.3| ) is illustrated in Fig. |2| for an arbitrary choice of tbbare^) = 560 MeV. It is seen 
that the effect of increasing the width is to change the slope of the real part of Tq (s), -Rq(s) 



near threshold from positive to negative, which contradicts experiment [44 . Note that the 
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FIG. 2. Predicted Real part of tttt 1=0 s-wave amplitude using regularization 
with different constant widths according to Eq. fl2.3|) . The widths are 50 MeV 
(solid), 100 MeV (dots), 200 MeV (dashes), 300 MeV (long-dashes), 500 MeV 
(dot-dashes). Here, weareC ") = 560 MeV. 

unitarity bound |i2g(s)| < \ is violated not too far away from threshold. Theoretically, it is 
most natural to use instead of an arbitrary constant, the "running" perturbative width, 



T(s) 



4m?. 



™bareI°"J 



(2.10) 



WttF^ V 

as was tried also in 0j. A plot of the real part of the resulting amplitude Rq(s) is shown 
in Fig. H and is seen to badly disagree with experiment. This is due to the large value 
of the perturbative width r [s — m^ ARE (a)}. ^ ^ s amusing that the somewhat arbitrary 
modification in the last factor of Eq. (|2.10| ): 



™bareW 



2 






s - ml 







(2.11) 



greatly improves the agreement near threshold, as is also shown in Fig. [3[ However somewhat 
beyond threshold the amplitude also starts to deviate greatly from experiment. Thus the 
prescription Eq. (|2.11|) does not completely solve the problem, but may help fitting to 
experiment if the effects of other possible particles are also taken into account. The effect 
of different values of h-bare (& ) in Eq. Q2.9| ) is illustrated in Fig. |]for this scheme. 



B. K-Matrix unitarization 



We can force unitarity at all s for the scalar partial wave amplitude, T®(s) by taking the 
tree amplitude TQ tree (s) given in Eqs. ( |2.8|) and ( |2.9|) to coincide with K(s) in Eq. (|i~l~l ). To 
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FIG. 3. Variations on SU(2) linear sigma model prediction for Real part of 7T7T 
1=0 s-wave amplitude. The dashed curve is the "current algebra" result Eq. 
(|2.2j ), the solid curve uses Eq. (|2.3| ) with the constant perturbative width calcu- 
lated from Eq. ( [2.10|) and the dotted curve uses the regularization prescription 
outlined around Eq. (|2.11|). Here, mBARE(o") = 560 GeV. 




0.25 0.3 0.35 0.4 0.45 0.5 



FIG. 4. Effect of different values of ^bare(c") on the generalized "running 
width" prescription, outlined around Eq. ( |2.11|) , on the SU(2) linear sigma 
model prediction for Real part of the tttt 1=0 s-wave scattering amplitude. We 
show the curves for (in MeV) it-bare^ ) = 300 (dots), 400 (solid), 500 (dashes), 
600 (long dashes) and 800 (dot-dashes). 
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see what is happening first consider Tq tree (s) to be small (for example near the irn scattering 
threshold). Then, in this single channel case, 

S° = 1 + 2zT ° = 1 + 2zT ° trce (s) + (2.12) 

so that Tq (s) starts out as To tree (s), which is a reasonable approximation. A presumably 
better approximation is obtained by including more terms in an expansion of the denomi- 
nator: 



T° r 2 

rrO( „\ _ 1 tree _ T 1 , - T , f-rpO \ , 

1 0\ b ) — ^ -rpo — J 0tree 1 + " tree + I 11 tree I ^ 
1 ~~ 11 tree 



(2.13) 



As observed in 0] this has the structure of a bubble sum in field theory. However, in the 
present case one is working with the partial wave, rather than the invariant, amplitude so 
there is no integration over intermediate state momenta. Of course, in either case, crossing 
symmetry is lost. While TQ tree is gotten from a crossing symmetric invariant amplitude, it 



is unlikely that the specially iterated amplitude Eq. (|2.13|) can be gotten in this way. The 
advantage of the method is that it guarantees unitarity. If, as is common, Tg tree (s) starts 
getting too large, Tq will be chopped down to size. For example if Tq tree gets very large: 

Sg— T °^z. (2.14) 

The real part R® of Tq vanishes in such a case while the imaginary part I® — > 1. In 
particular this occurs, as we see from Eq. Q2.8|), at the pole of TQ tree , where s = rri^ ARE . 
With the tree level amplitude of Eq. (|2.8|), the unitarized S-matrix takes the form: 



5 o/ x = [1 + ia(s)] [ml ARE (a) - s] + i/3(s) 
[1 - ia(s)} [m| ARE (a) - s] - i(3{s) ' 

where a(s) and /3(s) are given in Eq. ( |2.9| ). Eq. ( |2.15| ) is sufficient for comparing the 

predictions of the model [which contains the single unknown parameter ^bare^)] with 

experiment. However it is also of interest to rewrite the amplitude so that it looks more like 

a conventional resonance in the presence of a background. Manipulating Eq. ( p. 15 ) gives 

the factorized expression 

a, 



S (s) _ e 2i6 bg (s) rt (s)-S + t(3'( S ) 



where, 



tan [6 bg (s)] = a(s), 

m' 2 {s) = m BARE (a) + 



l + a 2 {s)' 



1 + a z (s) 
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This has the desired form although it should be noted that m! and ft' are both s-dependent. 
The T-amplitude which follows from Eq. ( |2.16| ) and Eq. ( |1 . 1| ) is the sum of a background 
term and a modified resonance term 



T°(s) = e iS »^ s hm6 bg + e 2 ^ (s) 9 ^ -. (2.18) 

m — s — ift' 



It is important to observe that the resonance mass and width (corresponding to a pole 
in the complex s plane) are shifted from their bare values. These new values should be 
obtained from the complex solution^, z a of: 

m' 2 (z) - z - ift' (z) = 0. (2.19) 

We may choose to identify the physical mass and width of the a from 

ml - im a T a = z a (2.20) 

One should keep in mind that the resonance term is no longer precisely of Breit-Wigner 
form. 

A plot of the real part R® (s) of Eq. ( 2 . 1 8|) is presented in Fig. [5] for the choices of 



(the single parameter in the model) wbare (c) = 0.5 GeV, 0.8 GeV and 1 GeV. It is seen, 
as already noted in [|]], that there is reasonable agreement with experiment up to about 
y/s = 0.8 GeV if m B ARE (c) lies in the 0.8 to 1 GeV range. Beyond y/s = 0.8 GeV, the 
effects of the /o(980), which does not appear in the two flavor model, are clearly important. 
Also, the unitarized curves for 7?ibare (c) in the 0.8 to 1 GeV range give a reasonable looking 
description of the threshold region, as opposed to the conventional unitarization scheme of 
Eq. (p|). 

Fig. P shows how the K-matrix unitarization works in detail by comparing R° (s) with 
Ro(s) tree . It is seen that Rv tree already violates the unitarity bound at ^Js = 0.43 GeV. 

Since we are regarding the K-matrix unitarization as a method of approximating all 
the higher order corrections to the nir scattering amplitude, it is clear that the quantities of 
physical significance should not be the bare mass and width of a but rather the pole mass m a 



ft For complex arguments the In function in a(z) is chosen to have an imaginary piece lying between 
— m and in. 

ttA different definition of the resonance mass and width was used in Q but the numerical results 
are close to each other. 
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m BARE (a) (GeV) 0.5 0.8 1 

r B ARE(cr) (GeV) 0.311 1.58 3.22 

m a (GeV) 0.421 0.458 0.449 

T a (GeV) 0.202 0.476 0.624 

z a GeV 2 ) 0.177 - 0.085 i 0.210 - 0.218 i 0.202 - 0.281 i 

a a (GeV 2 ) -0.015 + 0.078 i 0.088 + 0.169 i 0.158 + 0.188 i 

b a -0.420 + 0.443 i - 0.324 + 0.704 i -0.274 + 0.753 i 
TABLE I. Physical o parameters in the two flavor linear sigma model. 



and width r CT defined (with a usual convention) by Eqs. ( |2.19|j and fl2.2U| ). These quantities 



were obtained numerically and are given in Table 1 for the three choices of ttj-bare (c) used 
above. Evidently there are very substantial shifts of the bare mass and the bare width. The 
physical sigma pole mass is around 0.45 GeV while the pole width is around 0.5-0.6 GeV 
for mBARE (c) in the 0.8 to 1 GeV range. 

In the present model we may qualitatively understand the decrease in the a mass and 
also width by noting that a(s) and /3(s) vary slowly with s. If they are taken to be constant 
the physical mass m a would coincide with m' in Eq. ( |2.17|) and the physical quantity m a Y a 



with f3' in Eq. ( |2.17| ). Thus the negative sign of the mass shift arises since the background 
piece of the amplitude ot{s) is negative. A rather rough estimate may be made by evaluating 
a(s) and j3(s) for m n = and s small. Then one finds j3 — > 3m ?/ LT ^ a ^ > while a — > t-^—t\ - 

lb7rl< ir m BARE^^ 

It is interesting to note that our calculated amplitude Tq (s) can be reasonably well- 
approximated as 



T o (s) = + ba7 (2.21) 

S Zfj 



where the two complex numbers a and b are given in Table 1 for different choices of mBARE (c) . 
Since this simple pole dominated form reasonably fits experiment until the 700-800 MeV 
range it is not surprising that various determinations of m a and r CT in the literature are 
roughly similar to the ones in Table 1. Often the a parameters are stated in terms of z$. In 
the case where m a = 0.458 GeV we have 22 = 0.517 - z0.240 GeV. This may be compared, 
for example, with a treatment using a non-linear sigma model and including the p meson 
1% . That treatment gave a best fit for = 0.585 — i0.170 GeV. When it was refit (13 



without the p it yielded z? = 0.493 — i0.319 GeV, which is closer to the value in the present 
study (wherein, of course, spin 1 particles have not been included). 
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FIG. 5. Comparison with experiment of Real part of the I=J=0 txtx scattering 
amplitude in the SU(2) Linear Sigma Model, for mBARE(o") = 0.5 GeV (dots), 
^bare(o") = 0.8 GeV (dashes) and ^bare^) = 1 GeV (solid). Experimental 
data [44] are extracted from Alekseeva et al (squares) and Grayer et al (triangles). 
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FIG. 6. Comparison with experiment of Real part of K-matrix regularized 
I=J=0 7T7T scattering amplitude with the diverging real part of the tree ap- 
proximation, for me are (°") = 0.8 GeV. 
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III. SCATTERING IN THREE FLAVOR LINEAR SIGMA MODELS 



Here we study the pseudoscalar meson scattering amplitudes in the three flavor linear 
sigma models discussed in the Introduction. We shall restrict attention to the J=0 elastic 
scattering amplitudes of 

TTTT > 7T7T, 

ttK > 7tK, 

nr] — > nr], (3.1) 
which contain the scalars of the model - (a , k, a, a') - in the direct channel. 



Exactly the same K-matrix unitarization scheme of Eqs. (1J-) and ( |1.2|) which was used 
in the two flavor model will be employed. In particular, no special assumptions about the 
interplay of the a and a' resonances in tttt scattering will be made. The tree amplitude will 
simply be identified with K in Eq. ( |1 . 1| ) . The interesting question about the treatment 
of tttt scattering is whether it can fit the experimental data, given the complicated strong 
interferences between the a, a' and contact term contributions. The interesting question 
about the ttK scattering concerns the properties of the k meson in the present model. 
Finally the tttj scattering is of methodological interest. This is because the well-established 
a (980) resonance is expected to appear in a very clean way, lacking interference from a 
strong contact term (or even the possibility of potential interference when vector mesons are 
added to the model), as explained, for example in p4[ . 



We will first carry out the calculations using the standard renormalizable form of the 
three flavor linear sigma model. This is characterized by the potential in Eq. ( |1.21| ). Then 



the whole model is extremely predictive! After using as input the well-established masses of 
the pseudoscalar nonet and pion decay constant [Eq. (|1.14j )1 there is only one quantity left to 



choose in order to specify the scattering amplitudes. This one quantity may be taken to be 
the bare a mass, wt. B are(o")- The corresponding values of ttjbare^') and 6 S are given in Fig. 
[I]. We shall also carry out the calculations for the most general chiral symmetric potential. 
This allows ttj-bare^') and 9 S to be freely chosen, which is helpful for fitting experiment. As 
a possible justification for using a non-renormalizable potential we mention that the model 
is an effective one rather than the underlying QCD. (It may be considered, for example, to 
be a Wilson-type effective low energy Lagrangian. While non-renormalizable terms in the 
potential are technically irrelevant they play a part in establishing the spontaneously broken 
vacuum state and should be retained). In any event the extra parameters are being added 
in a chiral symmetric way. 
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A. 7T7T scattering 



The elastic amplitude for the three flavor linear sigma model in the tree approximation 
was given in Eq. ( |1.18| ) above. Calculating the I=J=0 partial wave amplitude as in section 
IIA gives a result which is a straightforward generalization of Eq. 



T 0tree(s) = COS^ 



a(s) + 



™BARe( ct ) - s 



+ sin 2 -?/' 



a{s) + 



™BAReK) - s 



(3.2) 



where a(s) and /3(s) are respectively gotten by replacing m B ARE( cr ) — ► ^bare^') in ®(s) 
and P(s) of Eq. fl2.9|) . The formula (|3.2|) evidently represents a sum of the o and a' related 
contributions, weighted by coefficients depending on the bare a — a' mixing angle ip. As 
before, we investigate the unitarized amplitude based on Eq. (|3.2|): 



™ = T^&Y (3 ' 3) 

which is being interpreted as an approximation to including the effects of the higher order 
corrections. 

First, consider the renormalizable model. To see how well it predicts the interesting 
I=J=0 amplitude we may simply plot the real part of Eq. (|3.3Q , Rq(s) against s for various 
choices of the single undetermined parameter ?71bare(o" ) • (Since we are working in an elastic, 
exactly unitary approximation, the imaginary part of the amplitude directly follows from 
the real part). Clearly the effect of the a' contribution is small when the mixing angle i/j 
is small. Then one returns to the two flavor case discussed in section II. Recalling that 
ip « 9 S + 54.7° and referring to the second of Fig. p] shows that this will be the case when 
^bare(o") is around 0.6 GeV. Rq(s) for this case is shown as the dotted line in Fig. [7[ The 
experimentally derived points are included for comparison. This looks similar to the 
curves for the two flavor case shown in Fig. [|, except that a sharp blip appears at about 1.09 
GeV corresponding to the a' pole in Tg tree . Note, as can be seen from Eq. ( |2.14| ), that the 
real part Rq{s) will vanish both at \/s = h-bare(c" ) and y/s = rn B ^ RE (a') . We saw in section 
IIB that the strong interference with the background appreciably changes the position of 
the physical pole in the complex s plane so, while useful for understanding Fig. |7|, these 
zeroes of R® do not give true parameterizations of the pole position. 

We have already learned that the choice ^bare^) = 0.6 GeV is too low for a good fit 
to R%(s) in the region of \fs up to about 0.6 GeV. Increasing 771bare(c" ) improves the fit 
to this region and also allows the effects of the a 1 to come into play. This is shown as the 
solid line in Fig. [7| which corresponds to the choice ?77.bare(0" ) = 0.73 GeV. Unfortunately 
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FIG. 7. Comparison with experiment of Real part of the I=J=0 tttt scattering 
amplitude in the renormalizable SU(3) Linear Sigma Model, for tit-bare (o~ ) = 
0.45 GeV (dashes), mBARE(c) = 0.6 GeV (dots) and ttibare(<7 ) = 0.73 GeV 
(solid). 

(as expected from Fig. |5]) this choice is still not high enough for a good fit in the low energy 
region. Furthermore, the structure which should correspond to the experimental /o(980) 
resonance has been pushed too high. Increasing /71bare(0") further will, as the first of Fig. 
|l] shows, push this structure even higher. One still must check to see if lowering tub are (c ) 
below 0.6 GeV can work. The dashed line in Fig. [?| shows Rq(s) for me are (°") = 0.45 GeV. 
In this case the structure near the /o(980) is in the right place but much too narrow. The 
structure at lower energies is also very badly distorted. 

By examining the evolution with increasing ttibare (c ) of the three curves in Fig. [7], one 
sees that it is not possible to get a good fit for any choice of tbbare^)- Nevertheless the 
qualitative prediction of Rq(s) is clearly a sensible one. It is therefore tempting to see if 
there is an easy way to fix up the fit. 

In a sense, the difficulty in obtaining a good fit arises because only one parameter - taken 
to be mBARE(c) - is available for adjustment to give agreement with a rather complicated 
experimental shape. The easiest way to proceed is to modify some parameters involved 
in the calculation. If a parameter to be varied is arbitrarily chosen there is however a 
danger of breaking the chiral symmetry relations intrinsic to the model. For example, 
suppose we choose to vary the coupling constant of the bare sigma to two pions. This three- 
point coupling constant, as mentioned in section I, is related to the masses of the particles 
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a a 1 

m B ARE ( GeV ) °- 847 L300 

r B ARE (GeV) 0.830 4.109 

m (GeV) 0.457 0.993 

r (GeV) 0.632 0.051 



z (GeV 2 ) 0.209 - i 0.289 0.986 - i 0.051 

a (GeV 2 ) 0.167 + i 0.210 0.053 - i 0.005 
b -0.248 +0.856 i 

TABLE II. Physical a and o' parameters obtained from best fit using non- 

renormalizable SU(3) linear sigma model. 

involved (see Appendix). Then, changing it without changing the masses will break the 
underlying chiral symmetry. Of course, we have written the formula for the tree ampliutde, 
Eq. ( |1.18p in such a way that the correct relations for the 4 and 3 point coupling constants 
are automatically taken into account for any choice of the contained parameters ^bare(c), 
^bare(<^ / ) and 9 S . In fact these three parameters may be freely chosen in the linear sigma 
model Eq. ( |1.5| ) with an arbitrary chiral invariant potential, Vq. It is only by restricting Vq 
to be renormalizable that one can relate weare^') and 9 S to w-bare (o* ) • Thus we can freely 
vary meARE^') and 9 S in addition to wibare("" ) if we choose to obtain the tree amplitude 
from the non- renormalizable model of Eq. flL5|) . 

In such a model it is easy to fit the experimental data for Rq{s). A best fit obtained 
using the MINUIT package is shown in Fig. |^. It corresponds to the parameter choices 
^bare(c) = 0.847 GeV, wt-bare(c"') = 1-30 GeV and ip = 48.6°. The physical masses and 
widths are obtained, as in Eq. ( [2.20|) in the two flavor case, from the pole positions in the 
complex s plane. These, together with the residues at the poles, are listed in Table [T|. 

For orientation we first note that the parameters describing the lower mass scalar, a are 
in the same range, as expected, as the parameters of Table | which give good fits to the low 
energy data using the a in the two flavor case. In fact the masses are very close to each 
other but the effect of the additional flavor requires a somewhat greater width parameter. 
The contribution of the <r-pole to Tq(s) is read off as: 

a a _ ctfj _ -0.167 — 10.210 

s - z a ~ ml- s- im a T a ~ 0.209 - s - i0.289' ' ' ' 

Note that this form is very different from a pure Breit-Wigner form which would require 
the numerator to be 0.289 GeV 2 . This illustrates (as does the two flavor case) the importance 
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FIG. 8. Comparison of our best fit for the Real part of the I=J=0 tttt scattering 
amplitude in the non-renormalizable SU(3) Linear Sigma Model with experi- 
ment. 

of the interplay between the resonance and the "background" . It also illustrates the possible 
difficulty of trying to get properies of the a from experiment with the use of a pure Breit 
Wigner approximation. Note (again) from Table [TT] that the physical mass and width, m a 
and IV respectively, are very much reduced from their bare values. The detailed mechanism 
is evidently similar to what we described in section II for the two flavor case. It seems 
reasonable to consider the physical values of m a and IV to be the ones which are significant. 

The well-established resonance /o(980) will be identified with the a'. The contribution 
of the a' pole to Tq(s) is read off as: 

a a < _ -a a > _ -0.053 - i0.005 

s - z a > ~ m 2 a ,- s- imjTj ~ 0.986 - s - z0.051 ' ^ ' ' 

In this case the form of a pure Breit-Wigner would require that the numerator be +0.051. 
To a reasonable approximation this holds except for an overall sign. Now reference to the 
formula Eq. ( |2.18j ) for a Breit Wigner with a background, shows that the background phase 
<5b g = | must be supplying this negative sign. Clearly the negative sign is required by the 
experimental data showing the real part -Rq(s) to be negative before and positive after the 
resonance. It was noted |12| that this is an example of the well-known Ramsauer-Townsend 



effect in scattering theory. It is also interesting to observe from Table |I| that the bare mass 
of the a 1 is substantially shifted down from 1.300 GeV (where a zero of Rq(s) remains, as 
previously discussed) to about 1 GeV. The bare width is even more substantially shifted 
from about 4 GeV to 50 MeV! 
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FIG. 9. Comparison of (left) Real and (right) Imaginary parts of pole approxi- 
mation Eq. ( |3.6| ) (dashed line) with those of our predicted amplitude Eq. (3.2) 
(solid line). 



One might wonder whether the simple pole dominance approximation Eq. ( |2.21| ) for the 
two flavor case can be generalized to this more complicated three flavor case containing two 
poles. It turns out to be true; the prediction of our model can be numerically approximated 
by the sum of the two pole terms and a suitably chosen constant: 

CL„ drrl 



T°(s) = R° (s)+iI (s 



+ 



+ 6, 



(3.6) 



where the numbers a a} a a >, z a , 



s - Z a S-Z a t 

z a i and b are listed in Table ||. This is illustrated in Fig. 



|9] for Rq(s) and Iq(s). In these figures Eq. ( |3.6|) is being compared with our prediction in 
Fig. |8|. Of course, there is no reason to use Eq. (|3T6| ) instead of the more accurate and 



complicated formula Eq. (37 3") but it nicely shows the dominating effect of the poles. The 
pole approximation is seen however not to be very accurate near threshold. 

Incidentally the deep dip in Iq(s) at what we have found to be the a' physical pole 
position also represents the Ramsauer-Townsend effect. This appears in -Rq(s) as a "flipped" 
resonance curve, as discussed above. Actually this Ramsauer-Townsend phenomenon can 
be pictured in an alternative manner. If we consider i?[} tree (s) corresponding to two "bare" 
resonances, one following the other, we see that there must be a point in between them where 
i?Q troe = 0. Then Eq. ( |3.3|) shows that, after K-matrix unitarization, R$ will also vanish 
at this point. This point appears visually as the pole position zero of a "flipped" standard 
resonance curve! In the Ramsauer-Townsend interpretation the flipping is interpreted as 
a background phase of |. Our explicit determination of the pole positions for the sigma 
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model amplitude shows that this is the pole which captures the dynamics of the /o(980). Its 
narrow width is seen to be the result of its getting "squeezed" between two nearby "bare" 
poles by the unitarization in this model. 



B. 7r K scattering 

We are interested in the 1=1/2, J=0 scattering amplitude in order to investigate the 
properties of the k resonance in the direct channel. The tree level amplitude involves k 
exchanges in the s and u channels, a and a' exchanges in the t channel as well as a four 
point contact term. The relevant tree level invariant amplitude may be written as: 

A 1/2 (stu) = -# (4) + 3 9 * K ^ 1 dlxir 9a^9aKK 9 'a' Vtt 9 'a> K K 

(3.7) 

where s, t and u are the usual Mandelstam variables. The four point contact interaction 
and the bare three point coupling constants shown are listed in the Appendix. As for the 
cases of a and a' we have put a subscript BARE on the k mass to indicate that the location 
of the physical pole after unitarization may come out different from this. The scalar partial 
wave tree amplitude is next defined by 



p(s) ^ dcos9A 1/2 (s,t,u). (3. 



T l/2 
1 Otree 

Note that p(s) was already defined by Eqs. ( |2.6|) and (|2.7|) . The specific formula for Eq. 
( |3.8|) in the present model is a bit lengthy and is shown in the Appendix. 

According to our plan we do not introduce any new parameters for unitarization and 
simply write 

T l/2 

rpl/2 _ -'Otree /o n\ 

~ 1 _ ,-T 1 / 2 ' 

1 11 Otree 

which is related to the corresponding S-matrix element by Eq. (|1 . 1|) . 

As mentioned in the Introduction the value of ?71bare(^) is independent of whether or 
not the chiral invariant potential in Eq. ( |1.5| ) is renormalizable, but depends only on the set 
of input parameters [e.g. Eq. (1.14)]. This may be seen from the equation 



"4™w = F «f- F ; m * (3.io) 

CK — fir 



which follows from Eqs. ( |1.9| ), ( 1.12|) and ( |1.17| ) in the isotopic spin invariant limit. This 



means that there are no new unknown quantities beyond those used in the fit to the tttt 
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scattering amplitude above. However we observe that the predicted value of /beare^) is 
very sensitive to the difference Fk — F w . Actually the choice of input parameters given in 
Eq. (|1.14Q results in a somewhat too high prediction for Fk, as mentioned before. 

One might therefore wonder whether the choice of input parameters in Eq. ( p..l4j ) un- 
fairly biases our treatment of ttK scattering by giving a too small value for wibare(^) ■ 
In order to check this we will also consider the slightly different choice of input param- 
eters 0(m,r, rriKi fTty, F n , Fk)- This will not affect the nit scattering results in the non- 
renormalizable model just discussed. We first choose Fk = 1.16^ which is slightly smaller 
than the physical value but has the advantage that it gives ?71bare(k) = 1.3 GeV which 
yields a zero for T 1/2 (s) at 1.3 GeV, in agreement with the experimental data. For this new 
input set we also have explicitly checked that there is still no possibility of getting a good 
fit to 7T7T scattering in the renormalizable model. 

With this choice of input and other coupling constants taken in agreement with those 
found in the best fit to tttt scattering we have the prediction for the real part of the amplitude 
R (s) shown in Fig. [It]. The experimental data j45], which start around 0.83 GeV and go 



to about 1.6 GeV are also shown in this graph. It is seen that the prediction from the linear 
sigma model agrees with the data from about 0.83 GeV to about 0.92 GeV. However at higher 
energies the predicted curve lies much too low until about 1.35 GeV and thereafter seems 
to completely miss the structure which is usually associated with the _K"q(1430) resonance. 

Fig. |Tt] also shows the predictions for the cases when wbare(^) = 11 GeV (correspond- 
ing to Fk taking its experimental value) and wt,bare(^) = 0.9 GeV (corresponding to the 
input choice of Eq. (|1.14j) . These are in worse agreement with experiment and also seem to 
miss the _K"q(1430) structure. 

As in the two flavor tttt case, which also contains only a single direct channel resonance 
we have found that the predicted amplitude is fairly well approximated as the sum of a pole 
term and a constant: 

r 1/2 ( s )^^^ + 6 K . (3.ii) 

s - z K 

The values for z K , a K and b K corresponding to the three different choices of input param- 
eters are shown in Table [TTJ. Again we identify the physical mass and width by 



ni 



l-im K T K = z K . (3.12) 



^We then predict m„ ~ 0.53 GeV rather than the experimental value of 0.547 GeV 
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FIG. 10. Comparison of our prediction for the Real part of the I=§, J=0 irK 
scattering amplitude in the non-renormalizable SU(3) Linear Sigma Model with 
experiment. The curves correspond to h-bare(^) = 1-3 GeV (solid), 1.1 GeV 
(dashed) and 0.9 GeV (dotted). The experimental data are extracted from [45]. 

It is notable that the pole position mass is always close to 800 MeV regardless of the choice 
of mBARE(ft)- Furthermore the widths obtained from Eq. ( |3.12| ) are substantially reduced 
from their "bare" (tree level) values, but are more sensitive to the choice of h-bare(^)- 

All in all, the properties of the k obtained here are very analogous to those of the a 
in either the two or three flavor treatments of titc scattering. Compare with Fig. [|, for 
example. 

It does seem that the pole mass, Eq. (|3.12| ), of the k is a good indication of the energy 
region where it provides a reasonable fit to the data. It also seems clear that the physics 
associated with the higher mass _K"q(1430) is not being taken into account in this model. 




C. irr/ scattering 



The tree level invariant amplitude takes the form 

1 



A 1 (s,t,u) 



+ 



-9 (4) + Q 2 

QanTiQar}!} 

^bare(o-) -t 



m 



+ 



bareM -8 m 

9 'a' irw 9 'a' r/rj 



2 

BARE 



(a ) - u 



™BARE<y) 



t 



(3.13) 



where the four point contact term — g^> as well as the three point coupling constants are 
listed in the Appendix. Other conventions are the same as above. Similarly the scalar partial 
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m B ARE(*0 ( GeV ) °' 9 L1 L3 

6 P 1.4,2.5° 1.23,-4.6° 1.16,-8.8° 

r B ARE(K) (GeV) 0.403 1.138 2.35 

m K (GeV) 0.799 0.818 0.798 

r K (GeV) 0.257 0.461 0.614 

z K (GeV 2 ) 0.639 - i 0.205 0.669 - i 0.378 0.637 - i 0.490 

a K (GeV 2 ) -0.043 + i 0.190 0.096 + i 0.340 0.263 + i 0.378 

b K -0.438 + i 0.420 -0.419 + i 0.660 -0.357 + i 0.800 

TABLE III. Physical k parameters obtained in the non-renormalizable SU(3) 
linear sigma model for different values of tj.bare( k ) which result from different 
choices of input parameters. 

wave amplitude is 

TLee = p{*) f\ '*xxB A 1 (s,t,u), (3.14) 
which is also listed in the Appendix. Again we unitarize by substituting this into the formula 

Tl = -j%g . (3.15) 

Since there is apparently no experimental phase shift analysis available for this channel, 
we will have to be content to just present our theoretical results and compare with the 
mass and width of the experimental ao(980) resonance. It was already noted that the 
renormalizable model [with the input Eq. (|1.14|) 1 yields the somewhat too low bare mass 
(which gets shifted down by unitarization) of 913 MeV. We will also present the results for 
the non-renormalizable model which gave a good picture of tttt scattering and for which 
we are still free to choose w^bare(^o)- A value w^bare(^o) = 1-100 GeV gives roughly the 
correct "physical mass" and the plot of the real part of Eq. ([3.15 ) for this choice is shown 



in Fig. |TT]. The result of the regularization is generally similar to the curves obtained for 
the a in tttt scattering and the k in ttK scattering. We have found in this case too that the 
predicted amplitude is reasonably well approximated by the sum of a pole and a constant: 



Tfo) = —*- + b ao . (3.16) 
The physical mass and width are found from 

m L - im ao r ao = z ao (3.17) 
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m BA RE(ao) (GeV) 0.913 1.100 

r B AR E (ao) (GeV) 0.129 0.381 

m ao (GeV) 0.890 1.013 

r ao (GeV) 0.109 0.241 

z ao (GeV 2 ) 0.793 - i 0.097 1.027 - i 0.244 

a ao (GeV 2 ) -0.065 + i 0.064 - 0.076 + i 0.200 

b ao -0.299 + i 0.204 -0.312 + i 0.408 

TABLE IV. Physical ao parameters in renormalizable (first column) and non- 
renormalizable (second column) SU(3) linear sigma model using corresponding 
best-fit parameters from tttt scattering. 

and the appropriate values for the two cases mentioned are listed in Table |V[ While Fig. |11 



seems to be just what one would expect for the real part of a resonance amplitude, Table [TV], 
as in the previous cases, reveals some interesting features. First, since a ao is clearly different 
from —lm(z ao ), the resonance is not a pure Breit Wigner. The location of the physical 
pole is close to the positive peak of Rq(s) rather than to its zero, as would hold for a Breit 
Wigner. Compared to the scalar resonance tttt and ttK channels we notice that there are 
smaller shifts going from 77Jbare(<^o) to m ao and from Tbare^o) to r ao . This is reasonably 
interpreted as due to less effect of interference with the background. This is manifest in 



the non-linear sigma model approach to tct] scattering |?4J] and can thus be understood as 
a consequence of the similarity of the non-linear and linear chiral models. In addition, we 
note that r ao is predicted to be somewhat larger than the experimental value [4~5] of 50-100 
MeV. Nevertheless, the prediction is qualitatively reasonable. 



IV. SUMMARY AND DISCUSSION 



We have treated the three flavor linear sigma model as a "toy model" for examining the 
role of possible light scalar mesons in the nn, ttK and 7rr) scattering channels. This is a highly 
predictive model which contains only one free parameter, which may to taken as ^bare(c), 
in the renormalizable case. If we give up renormalizability for this effective Lagrangian 
but maintain chiral symmetry in a straightforward way, t^bare^'X the sca l ar mixing angle 
9 S and ^bare(oo) may also be freely chosen, which is helpful for fitting experiment in 
the desired energy range of threshold to the 1+ GeV region. Our approach just involves 
computing the tree amplitude for each channel and unitarizing by a simple "K-matrix" 
prescription which does not itself introduce any new parameters. In general the unitarization 
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FIG. 11. Prediction for the Real part of the 1=1, J=0 tit] partial wave scatter- 
ing amplitude Eq. ( |3.15| ) in the non-renormalizable SU(3) Linear Sigma Model 
(parameters as in second column of Table |TV| ) . 
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has very important effects converting "bare" scalar meson masses and widths into "physical" 
ones. It turns out that there is not too much "wiggle room" in this procedure so that what 
results is characteristic of the model (and the unitarization scheme). This tightness comes 
from the demand that the starting tree amplitudes satisfy chiral symmetry restrictions. This 
means, as discussed in section I, that the four point contact interaction vertices are related 
to the three point interactions which are related to the particle masses (two point objects) 
which are related in turn to the one point terms (pseudoscalar meson decay constants). 
We chose the inputs to be the four pseudoscalar masses and the pion decay constant [Eq. 
( |1.14j )1. However the pseudoscalar mixing angle and kaon decay constant were not perfectly 
predicted so there is already a source of error present before even going to the scattering 
amplitudes. Nevertheless we investigated this point by considering an alternative input set 
obtained by using F% instead of m v and found that there was not much qualitative change 
for the scattering predictions. 

Our point of view in this paper is to see what are the results of computing in a relatively 
simple and natural model for the purpose of comparison with other (and possibly future) 
more elaborate treatments. It seems to us that the results are interesting and instructive. 
In the simpler two flavor case, which was applied in 0] to a lower energy treatment of tttt 
scattering, the results were already reasonable. Here we have, in section II, reviewed the 
two flavor case in a slightly different way as preparation for the more complicated three 
flavor case. We have also made some new comments and suggested an alternative "naive" 
unitarization procedure which might be handy for future studies. 

Table [V| contains a brief summary of the physical masses and widths of the scalar mesons 
predicted in the present model and discussed in some detail in section III. In the cases of 



the /o(980) and a (980) resonances comparison is being made with experimental values [06 



In the cases of the a and the k, which are less well-established experimentally, we have 
compared with the earlier computations of the Syracuse group ||l2|]l4l , |l5| which were based 
on a non-linear chiral effective Lagrangian treatment, including vector mesons. Many other 
authors ]3(J were led to similar predictions for the a while a similar prediction for the k was 
made by |p"T| . 

The predicted properties of the a and fo in the present model come from their role in titt 
scattering as discussed in section III A. It was found that the single parameter describing 
the renormalizable model could not be adjusted to give a reasonable fit to the experimental 
data. This could be done when the renormalizability condition was relaxed. Neither the 
physical a nor the physical fo are described by simple Breit Wigner terms. Both have masses 
and widths greatly reduced from their "bare" values by the unitarization procedure. The 
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a f k a 

Present Model 

mass (MeV), width (MeV) 457, 632 993, 51 800, 260-610 890-1010, 110-240 
Comparison 

mass (MeV), width (MeV) 560, 370 980±10, 40-100 900, 275 985, 50-100 

TABLE V. Predicted "physical" masses and widths in MeV of the nonet of scalar 
mesons contrasted with suitable (as discussed in the text) comparison values. 



light, broad a is somewhat lighter and broader than the comparison one obtained in the non- 
linear model |T2|] . (However when the vector meson contribution in the non- linear model 
was consistently eliminated [0 the a in that model also became broader and lighter). The 
/o obtained approximately looked like a Breit-Wigner in a background which has a phase 
5b g = | . This is known as the Ramsauer-Townsend effect in scattering theory. The fact that 
it emerges in the present model was noted to be explicable in terms of the region between 
two neighboring "bare" resonances getting squeezed by unitarization. 

The entries in Table [V] for the k mass and width require some explanation. The bare 
k mass and width in this model are uniquely predicted once the input parameters are 
specified, regardless of whether or not the potential is taken to be renormalizable. However 
the predictions of the k parameters are very sensitive to Fk (which measures the deviation 
of the vacuum from exact SU(3) flavor symmetry in this model). Thus we allowed different 
input sets yielding different bare k masses, as discussed in section IIIB. Whatever reasonable 
choice was made, the unitarization always brought the physical k mass down to around 800 
MeV. However the physical width is more dependent on this choice. Furthermore, as shown 
in Fig. [H], the k resonance can only explain the lower energy 7rK scattering data. This 
would be the analog of the SU(2) treatment of tttt scattering, where the a alone can provide 
a reasonable description of the low energy region. The k cannot explain the data in the 
region of the ^(1430) scalar resonance. In other words, we cannot explain the ^(1430) 
as the strange scalar of the usual linear sigma model treated with K-matrix unitarization. 

In the case of the ttt] channel there does not appear to be any experimental phase shift 
data, so we compare with experimental determinations of the ao(980) mass and width. 
The lower physical mass entry for the ao in Table [V] corresponds to the bare mass of the 
renormalizable model. It is somewhat too low but not very far off. This can be easily 
adjusted by using the non-renormalizable potential. The predicted width is somewhat too 
large but qualitatively reasonable. Clearly, the ao of the present model is describing the low 
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energy part of nr) scattering and should correspond to the ao(980) rather than the ao(1450). 

All in all, the three flavor linear sigma model with a general (non-renormalizable) chiral 
invariant potential and regularized by the simple K-matrix procedure can approximately 
describe the complicated nir scalar scattering amplitude as well as the low energy part of 
the Kn scalar amplitude and the a (980) nr] resonance. The _K"q(1430) and ao(1450) are 
"outsiders" in this picture and would have to be put in by hand to realize the higher mass 
scalar resonances in nK and irr] scattering. The picture is qualitatively similar to that 
obtained in treatments using the non-linear sigma model for ttti |12|], ttK [[Li] and m] [24 



scattering. The ao(980) and /b(980) seem to belong to the same multiplet as the controversial 
light a and light k. Of course, it is possible for particles with the same quantum numbers 
belonging to other multiplets to mix with them. 

There are several straightforward, but lengthy to carefully implement, ways to improve 
this treatment. Modified kinetic terms, as mentioned in section I, can be included to im- 
proved the fit to pseudoscalar masses and decay constants. Vector and axial vector mesons 
can be added to introduce more of the low-lying physical resonances which are expected to 
be important in the tctt and irK channels. Finally, inelastic channels can be included. 



V. SPECULATION ON SCALAR MESON'S QUARK STRUCTURE 

Up to this point we have reported the results of a straightforward and highly predictive 
treatment of the three flavor linear sigma model. Our original reason for pursuing this 
investigation was to check the results obtained in our treatment of meson scattering in the 
non-linear sigma model which contained additional particles and channels. That treatment 
used a different unitarization procedure in which crossing symmetry and unitarity were 
both approximately satisfied. (Actually in the study of direct channel scalar resonances, 
the crossed scalar exchanges are relatively small). We already noted that the locations and 
widths of the physical scalar states obtained in the linear model were qualitatively similar 
to those obtained in the non-linear model. Since the a, fo, ao and k all come out less than 
or about 1 GeV, and the scattering regions near the ao(1450) and the -fTo(1430) apparently 
must be described by fields other than those contained in the matrix M, the well known 
puzzle of the quark structure of these scalars comes to the surface. 

In this section we will make some speculative remarks on this controversial subject and 
introduce another toy model which may illuminate some of the issues. The puzzle, of course, 
is why, if the scalars are u qq states", they are considerably lighter than the other p-wave 
states and why the isovector a (980) is tied for being the heaviest, rather than the lightest, 
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member of the multiplet. 

Actually there is a lot of ambiguity in stating what the quark structure of a physical 
hadron means. Generally people think of the question in the context of a potential-type 
model wherein, for example, the p meson is made of a "constituent" quark of mass about 
300 MeV and a constituent anti-quark of the same mass. The idea is that the fundamental 
"current quarks of QCD" (with masses about 10 MeV) interact strongly with each other 
and with gluons to make the relatively weakly interacting constituents whose combined 
masses roughly approximate the physical hadron masses. Thus the quark structure really 
depends on the model used to treat the hadrons. At the field theory level of "current quarks" 
there is always some probability for extra qq pairs or other structures to be present. In the 
SU (3) i x SU (3)r chiral effective Lagrangian treatments, the quark substructure of the fields 
being used does not enter the formulation in a unique way. An infinite number of different 
quark substructures will give rise to the same SU(3)l x SU(3)r transformation properties 
for the mesons. This is apparent for the non-linear chiral model in which scalars are added 
to the pseudoscalar meson Lagrangian as "matter fields" in the usual manner |[43| . Then 
it is known that only the SU(3) flavor transformation properties of the scalars are relevant. 
However we found in our earlier study |l5j that the value of the scalar mixing angle suggested 



indirectly that the light scalars do have an important four quark component. Considering 



the properties of the heavier scalars a (1450) and _K"q(1430) suggested [25| that these states 
did not belong to a "pure" qq multiplet but to one which mixed with the lighter scalar 
multiplet. 

When it comes to the linear sigma model where the chiral transformations of the scalars 
are linked with those of the pseudoscalars in a natural way, there seems to be a feeling that 
the matrix M should describe a qq field. In fact, there are still an infinite number of quark 
substructures which transform in the same manner under SU(3)l X SU(3)r. It may be 
worthwhile to illustrate this for the specific cases of interest in the literature. 

The schematic structure for the matrix M(x) realizing a qq composite in terms of quark 
fields q a A(x) can be written 

M^ b = (q bA ) j l,^qaA, (5.1) 

where a and A are respectively flavor and color indices. Our convention for matrix notation 
is M^ b -> M$. Then AfW transforms under chiral SU(3) L x SU(3) R as 

M (1) -> U L M w U j R (5.2) 

where Ul and Ur are unitary, unimodular matrices associated with the transformations on 
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the left handed (gx = \ (1 + 75) <?) and right handed = \ (1 — 75) g) quark projections. 
For the discrete transformations charge congugation C and parity P one verifies 



C : M«^M« T , P: M«(x) -> #t(-x). (5.3) 
Finally, the {7(1)a transformation acts as q a L — » e w q a Li QaR ~^ e ~ zu QaR and results in: 

M {1) -> e 2 ^M {1) . (5.4) 



One interesting model 47] for explaining the scalar meson puzzle (at least insofar as 
the ao(980) and /o(980) states are concerned) is to postulate that the light scalars are 
"molecules" made out of two pseudoscalar mesons. The chiral realization of this picture 
would result in the following schematic structure: 

M^ b = e acd e^(M^ e (M^Y. (5.5) 



One can verify that transforms exactly in the same way as M^ l > under SU(3)l X 

SU(3)r, C and P. Under U{1)a it transforms as 

M (2) _ >e -M»' M (2) j (5. 6 ) 

which differs from Eq. (|5.4j) . 

Another interesting approach ^| to explaining the light scalar mesons was formulated by 
Jaffe in the framework of the MIT bag model. It was observed that the spin-spin (hyperfine) 
piece of the one gluon exchange interaction between quarks gives an exceptionally strong 
binding to an s-wave qqqq scalar state. Furthermore, this model naturally predicts an 
"inverted" mass spectrum of the type summarized in Table |V[ The scalar states of this type 
may be formally written as bound states of a "dual quark" and "dual antiquark" . There are 
two possibilities if the dual antiquark is required to belong to a 3 representation of flavor 
SU(3). In the first case it belongs to a 3 of color and is a spin singlet. This has the schematic 
chiral realization, 

l9 e = e^e EAB q T aA C^ l -±^q hB , 

r9 e = e^e EAB q T aA C- ll -^ qhB , (5.7) 

where C is the charge conjugation matrix of the Dirac theory. A suitable form for the M 
matrix is: 
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M^ 3 ) can be seen to transform in the same way as under SU(3)l x SU(3)r, C, P and 
U(1)a- In the second case the dual antiquark belongs to a 6 representation of color and has 
spin 1. It has the corresponding schematic chiral realization: 



T 9 -J 9 - c9 ah a T C~ l a 1 + 75 ru. n 



2 

R %,AB = R %,BA = e9<1 ' ^Ia^ L > ^~QbB, (5.9) 



where a^ u = ^ [7^, This choice leads to an M matrix 

M^ )f = (L% AB fR^ AB , (5.10) 

where the dagger operation includes a factor (_i) ( ^ 4 + 51 ' 4 . also transforms like 

and under all of SU(3)l x SU(3)r, C, P and U(1)a- The specific form favored by 
the MIT bag model calculation actually corresponds to a particular linear combination of 
and M^. Furthermore one can verify that in Eq. (|5.5| ) is related by a Fierz 

transformation to a linear combination of and M^. Thus only two of M^ 2 \ and 
are linearly independent. 
What is the significance of these remarks for construction of the general effective chiral 
Lagrangian used in this paper [Eq. ( [L.5|) 1? All that is required for M is that it transform 
like under SU(3)l x SU(3)r, C and P and that it carry a non-zero U(1)a "charge" 
which gets broken by the potential. The specific U(1)a transformation property does differ 
between the two quark realization and the four quark realizations (M^ 2 \ and 
M^) but this would just be absorbed, in the present work, by a different value for the 
parameter V4. Thus, if one knew nothing else about hadronic physics than the present toy 
Lagrangian, one would not be able to a priori easily discriminate among the possibilities 
M^ l > — , or in fact any others, for the underlying quark substructure of the scalars (and 
pseudoscalars) . Nevertheless, one might glance at the obtained scalar masses in Table 
and notice that there is an inverted physical mass spectrum. One might then decide to 
make a judgement on the constituent quark substructure by fitting the scalar spectrum to 
an Okubo type mass formula JOJ. This was done recently, for example, in section II of fl5| 
and suggests that the scalars are behaving roughly as composites of four constituent quarks. 
Roughly, this amounts to simply counting the number of strange constituent pieces in each 
state; in the four quark picture both / (980) and ao(980) have two. The combined effects 
of spontaneous chiral symmetry breaking and unitarization (presumably taking radiative 
corrections into account) appears to split the constituent structures of the scalars from the 
pseudoscalars, regardless of which current quark structure (i.e. choice of M) we start with. 
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However the true situation is likely to be more complicated. The present model does 
not appear to accomodate the a (1450) and 1^(1430) scalars as states belonging to M. 
These states would seem at first sight to be reasonable candidates for a nonet of ordinary 
qq scalars. Still it is a little puzzling that _K"q(1430) is not heavier than ao(1450). There are 
some other puzzles too but all can be qualitatively explained |25] if a qq scalar nonet mixes 
with a qqqq scalar nonet. If we want to realize such a scheme in the linear model framework 
it would be natural to introduce a Lagrangian with two different M matrices. Such a model 
seems to yield a variety of interesting dynamical possibilities which may lead to new insights 
and approximation schemes for low energy QCD. Thus it may be worthwhile to give a brief 
discussion here. 

Let us start with the field M^ x > which we shall simply designate M. At the kinematical 
level it represents a current-type quark antiquark operator. This is modified for both the 
pseudoscalar and scalar states by the (almost) spontaneous breakdown of chiral symmetry. 
For the scalars (which occur as poles in the physical region) there is an additional modifica- 
tion due to the unitarization required. Of course, the choice of the free parameters gives an 
"experimental" input to this process. The resulting scalars seem to be roughly consistent 
with a qqqq constituent-quark structure, as just discussed. Now consider adding a current 
type four quark operator which may be any combination of M^ 2 ' , or M^> (we could 

not tell the difference in an effective Lagrangian framework) and denote it by M' . Allow M' 
to mix with M. What happens? 

The Lagrangian which directly generalizes Eq. ( |1.5[ ) is written as 

C = -Itt (d,Md,M^) - (^M'^M't) - V (M, M') - V SB , (5.11) 

where Vo(M, M') stands for a general polynomial made from SU(3)l x SU(3)r [but not 
?7(1)a] invariants formed out of M and M'. Furthermore Vsb is taken to be the same as 
Eq. (1.7) since it is Tr + which "mocks up" the quark mass terms. Other physical 
particles (including glueballs) could be added for more realism, but Eq. (|5.11|) is already 
quite complicated. 

To get an indication of what kinds of questions might be answered, let us consider a very 
simplified approximation in which the quark mass effective term, Vsb is absent and where 
Vo is simply given by: 

V = -c 2 Tr (MM 1 ) + c 4 Tr (MM^MM ] ) + d 2 Tr (m' M' ] ) + eTr (hIM' ] + M'M f ) . (5.12) 

Here C4 and d 2 are positive real constants. The M matrix field is chosen to have a 
wrong sign mass term so that there will be spontaneous breakdown of chiral symmetry. A 
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pseudoscalar octet will thus be massless. On the other hand, the matrix field M' is being set 
up to have trivial dynamics except for its mixing term with M. The mixing is controlled by 
the parameter e and the e-term is the only one which violates U(1)a symmetry. Its origin is 
presumably due to instanton effects at the fundamental QCD level. (Other U(l) ^-violating 
terms which contribute to rj mass etc. are not being included for simplicity). Using the 
notations M = S + i<f) and M' = S' + i(f)' we may expect vacuum values 



S. 



The minimization condition 



aS° a , {S. 
= leads to 



lb 



(5.13) 



while (|| 



yields 



a" 



(3 = - T a 



— (c - 

2c 4 I 2 d 2/ 



(5.14) 



(5.15) 



In the absence of mixing the "four-quark" condensate (3 vanishes while the usual two 
quark condensate a remains. 

The mass spectrum resulting from Eq. Q5.12 ) has two scalar octets and two pseudoscalar 
octets, each with an associated SU(3) singlet. Each octet has eight degenerate members 
since the quark mass terms have been turned off. Let us focus on the 1=1, positively charged 
particles for definiteness and define: 



7T 



7T 



/+ 



q2 



12 



(5.16) 



Then the 2x2 squared mass matrix of n and it' is: 

2 



This has eigenstates 



t- e 

d.2 

e d 2 



(5.17) 



with masses 



m 2 (vr p ) = 0, 



m 



BARE 



(5.18) 



(5.19) 
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We put the subscript "BARE" on m 2 (j^'^j to indicate that it may receive non-negligible cor- 
rections from K-matrix unitarization as in our detailed treatment of the M only Lagrangian 
in the above. A possible experimental candidate for such a particle is the 7r(1300). 
Computing the axial vector current by Noether's theorem yields 



F 7 



2a\ 



' d 2 



(5.20) 



where a is given in Eq. ( |5.15 ). 

Notice that a term like d^n'^ does not appear in our semi-classical approximation. 
The 2x2 squared mass matrix of the scalars a and a' is 



4c 2 + 
2e 



6e 2 
d 2 



2e 
2do 



(5.21) 



The eigenstates are defined by 




cosu; -sinw 
sincj coso; 



(5.22) 



with 



tan2u; 



2d 2 - 4c 2 - Q 4- 



(5.23) 



The corresponding masses are 

(a p , a^j = 2c 2 



BARE 



d 2 



3e e 

d 2 H — — =F - (3sin2o; + csc2c<j) 
d 2 2 



(5.24) 



where the upper (lower) sign stands for a p , (a' p ^ . 

It is interesting to examine the masses of the degenerate octets in a little more detail. For 
orientation, first consider the case when the mixing parameter e vanishes. The usual "gg" 
pseudoscalars ir p are zero mass Goldstone bosons in this approximation. If 4c 2 > 2d 2 , a p , the 
original scalar partner of ir p lies higher than the degenerate u qqqq" scalar and pseudoscalar 
a' p and tt' p . When the mixing is turned on, a four quark condensate develops and the mass 
ordering is 



^BARE(a p ) > m B ARE(vTp) > m BARE (a' p ) > m B ARE(vr p ) = 0. 



(5.25) 



This is graphed, as a function of e, in Fig. [12] (with parameter choices c 2 = 0.25 GeV 2 , 
d 2 = 0.32 GeV 2 . In such a scenario, the qqqq scalar would be the next lightest after the qq 
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FIG. 12. Plots of mBARE (op) (solid), wibare («p) (dashed) and mBARE \x' p 
versus the mixing parameter e for the choice C2 = 0.25 GeV 2 and (L2 = 0.32 
GeV 2 . The highest lying curve is mainly a u qq" scalar, while the lowest lying 
curve is mainly a "qqqq" scalar. The excited pseudoscalar curve is in the middle. 



Goldstone boson. Each particle would be a mixture of qq and qqqq to some extent. For the 
given parameters the mixing angle remains small however because the denominator of Eq. 
( |5.23| ) is always negative and increases in magnitude as e 2 increases. Note especially, that 
due to the spontaneous breakdown of chiral symmetry, there is no guarantee that the lowest 
lying scalar is of qq type. Also note that it' is expected to be more massive than a' p . 

On the other hand, if the QCD dynamics underlying the effective Lagrangian is such 
that 2g? 2 > 4c2 we will get a mass ordering jit-bare^) > ^bare(tt') > ^bare ( a p) in which 
the four quark scalar appears heaviest. However, in this case we will definitely get a large 
mixing as e increases since the denominator of Eq. ( |5.23 ) starts out positive when e = 
and will go to zero as e is increased. Thus the next-to-lowest lying a p can be expected to 
have a large qqqq admixture. 

All of these remarks pertain to the meson current-quark type operators in the toy model. 
The important effects of unitarization (i.e. ttibare ~~ > m ) ar e likely, as in our earlier treat- 
ment, to favor an interpretation of the low lying physical scalars as being of four constituent 
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quark type in either case. 

The main lesson from our preliminary treatment of a chiral model with mixing is perhaps 
that even though the M fields carry "chiral indices" it is not easy to assign an unambiguous 
quark substructure. On the other hand there is a great potentiality for learning more 
about non-perturbative QCD from further study of the light scalars. Such features as scalar 
mixing (including the possibility of mixing with glueballs for the 1=0 states), four quark 
condensates and excited pseudoscalars may eventually get correlated with each other and 
with the experimental data on the scattering of light pseudoscalars. 
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APPENDIX A: COUPLING CONSTANTS AND PARTIAL WAVE AMPLITUDES 



For the scattering processes under consideration we will need the four-point pseudoscalar 
contact interactions and the trilinear scalar-pseudoscalar-pseudoscalar interactions. In iso- 
topic spin notation the relevant pieces of the Lagrangian are respectively: 

- £« = ^(tt • tt) 2 + ±g$KKn ■ tt + \g^ VW ■ tt + ... (Al) 

and 

- £s<t><t> = ^7=f {KT ■ TZK + /i.C.) + ^p<77T • 7T + g aKK (jKK + ^^a'lZ ■ 7T (A2) 

+ g a 'KK<y'KK + Kt ■ & K + g KKr) (kKt) + h.c.) + g KKri > {RKrf + h.c.) 
+ fi'ao^ao • ni] + gaol's*) ' + -^^VV + 9a m '^VV H 2 ar> 11 



i fl'c'njj i . ii, g<r"n"n' i i i / A o\ 

+ ~Y ±a VV + 9v'm'° VV + V V • ( A3 ) 

The trilinear couplings which do not involve three isoscalars are predicted in terms of the 
masses. These are given in |]37| and we present them here for completeness: 

9*K* = (™bare(«0 - m l) , 9*k v = ( cos6 p + 2V2sin0 p ) (m 2 - m 2 ARE («;)) , (A4) 

g K Krf = ^ F (2V2cos6 p - sin# p ) (m| ARE (/«) - m 2 ,) g aoKK = — (m| ARE (a ) - m 2 K ) , 
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(rn 2 BARE {a ) - m 2 ) , #„ ^ = (^IareM - mj,) , 

/ 2 /x 2\ \/2, / 



BARE 



a - to. 



#<rK*: = (cos# s + 2V2sin0 a ) (m|- - m| ARE (o-)) , 

S-a'^K = ^/^T (2v^cos0, - sin6> s ) (m| ; 



V) - . 



The trilinear coupling constants involving three isoscalars may depend on V4. For tti] elastic 
scattering we will also need: 

6, 



9. 



V2 



X-b s Y, g, 



V2 



X + a s Y, 



(AS) 



where 



x= (*, A 



2o>bare(^) + 2b 2 s m 2 BARE (a') - m 2 n - a 2 p m 2 v - b 2 p m 2 f - 12{2F K - F n )V 4 ] (A6) 
2 



- F w 



■V2a s b s ( 



m BAREl°"J — m BARE 



>')) - 12F^ 4 



48 , T T 



\V2J F w 



-V2a s b s (ml ARE ((r) - ml ARE (a')) - 24F W V 4 

l2 2 / \ 1 2 2 / A l2 2 22 
^BAReW + fl s m BAREl ff ) ~ b P m v ~ a P m v > 



p 2F K -F„ 

In these equations we have used the convenient abbreviations 



(A7) 



cos^p — y^sin^p 



v^2cos^ p + sin^p 



(A8) 



V3 ' p Vs 

with analogous expressions for a s = cosip and b s = simp in terms of 6 S . The contact coupling 
constants are then given by: 

4 



9n ] = ~p2 ( a s m BARE( ff ) + ^^BARe(^) ~ m l) , 

7T 



(A9) 



(4) 



y»7 



1 



2 / \ 2 2,22 / \ I 7 2 2 / /\ 

™BARe( k ) - m K ~ m n + a S m BAREl^) + & S ™BARe(c ) 



V2a s b s (m 2 BARE {a) - m BARE (a')) 

ctg bg ^2/2 / \ i 

~^9m m + -^Qa'rtr) + -pT a p ( m BAREl a o) — m i{ 



(A10) 
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Finally, the tree-level partial wave amplitudes for 7r.fr and irr) scattering are: 
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where q(s) and p(s) for each case are given by Eqs. ( |2.7|) and ( |2.6|) respectively. Furthermore 
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